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The Kardar-Parisi-Zhang universality class of stochastic surface growth is studied by exact field- 
theoretic methods. From previous numerical results, a few qualitative assumptions are inferred. 
In particular, height correlations should satisfy an operator product expansion and, unlike the 
correlations in a turbulent fluid, exhibit no multiscaling. These properties impose a quantization 
condition on the roughness exponent \ and the dynamic exponent z. Hence the exact values \ — 
2/5, z = 8/5 for two-dimensional and \ = 2/7, z = 12/7 for three-dimensional surfaces are derived. 
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Strongly driven dynamic systems offer some of the 
most intriguing realizations of statistical scale invariance. 
Hydrodynamic turbulence jlj or the growth of rough sur- 
faces ||] are two classic examples, which turn out to be 
deeply connected from a theoretical point of view. In 
such systems, a stochastic force r](r,t) generates long- 
ranged correlations of a fluctuating dynamic field - the 
local velocity v(r, t) of a fluid or the height h(r,t) of 
a surface. As typical differences h(ri,i) — h(r%,t) or 
v(ri,t) — v(r 2 ,t) increase with the spatial separation 
\t\ — r 2 |, the scaling properties of these fields are gener- 
ically more complex than those at a standard critical 
point. Indeed, the theoretical understanding of these uni- 
versality classes far from equilibrium is still fragmentary. 

The subject of this Letter is the simplest nonlinear 
model of stochastic surface growth, the famous Kardar- 
Parisi-Zhang equation 
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for a rf-dimensional surface q|. The driving term r)(r, t), 
which describes the random adsorption of molecules onto 
the surface, is taken to be Gauss distributed with corre- 
lations only over microscopic distances, 



r)(r, t)r)(r', t') = a 2 S(t - t')5(r - r') 
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The relation of this model to the theory of turbulence 
is manifest: Eq. (|l|) is formally equivalent to Burgers' 
equation 
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for the driven dynamics of a vortex-free velocity field 
v(r,t) = V/i(r,t) (with A = 1) (§. In a fluid, however, 
the driving force is correlated over macroscopic spatial 
distances. This leads to important differences in the scal- 
ing behavior Q , which are discussed below. 

A surface growing from a flat initial state h(r, 0) = 
develops height correlations with an increasing correla- 
tion length £ t ~ t 1 / 2 , which defines the dynamic expo- 
nent z i. A self-similar growth pattern, characterized 
e.g. by the height difference moments 



(with ri2 = ri — r2), emerges on mesoscopic scales 
a <C IrigJ <§C £t- For |ri 2 | ^o, the dissipation term v V 2 h 
in Eq. (|l) breaks the asymptotic scale invariance. In the 
scaling regime (^) , the height difference moments become 
stationary, i.e., independent of the correlation length £ t . 
They are characterized by a single critical index x > 0j 
the roughness exponent of the surface. The scaling re- 
lation x + z = 2 follows from the Galilei invariance of 
Eq- (0) @ • For d = 1, one can show Eq. (|J) to be valid 
with the roughness exponent x — 1/2, equal to that of the 
linear theory (A = 0) In higher dimensions, however, 
little is known analytically. For d > 2, the rough state of 
the surface exists only if the rescaled driving amplitude 
Aq = (T 2 A 2 /^ 3 exceeds a finite threshold value X 2 Less 
rigorous theoretical arguments predict an upper critical 
dimension d > < 4 beyond which Kardar-Parisi-Zhang 
surfaces are only logarithmically rough (x — 0) even in 
the strong-coupling regime X 2 > X 2 ||. The numerical 
results presently available are consistent with Eq. (Q). 
Extensive simulations yield x ~ 0-39 for d = 2, x ~ 0-31 
for d = 3, and smaller positive values in higher dimen- 
sions, which are less reliable 1 10 1 . 

It has remained a challenge for theorists to calculate 
the rough asymptotic state of Kardar-Parisi-Zhang sur- 
faces for d > 1 exactly or in a controlled approxima- 
tion. In particular, standard perturbative renormaliza- 
tion about the linear theory fails to produce a fixed point 
belonging to this regime (ll) - a notorious difficulty fa- 
miliar from the theory of turbulence. In this Letter, a 
quite different approach is taken. Guided by numerical 
and experimental results, I make a few qualitative as- 
sumptions, namely the existence of an operator product 
expansion (||) and of a stationary state ([ll]) that is di- 
rected (i.e., has no up-down symmetry). These assump- 
tions turn out to constrain severely the possible solutions 
of Eq. (Q). In particular, they naturally lead to a quan- 
tization condition for the roughness exponent: 
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where ko is an odd integer for d > 2. Comparing with 
the above numerical estimates 13 and using the relation 
X + z = 2 then gives the main result of this Letter: the 
exact values x = 2/5, z — 8/5 for d = 2 and x — 2/7, z 
12 7 ford = 3. 

The fundamental observables describing the equal-time 
surface configurations are the (connected) correlations 
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(the dots denoting the disconnected parts). The height 
probability distribution P t ({h}) obeys the functional 
Fokker-Planck equation 
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where J(r) = vV 2 h{r) + (A/2)(V/i) 2 (r) is the determin- 
istic part of the current. 

In the scaling regime (a <C \rjj\ <C £ t for i,j = 
1, . . . , n), the correlation functions (p|) will generically be- 
come singular as some of the points approach each other. 
For d < d>, these singularities are assumed to follow 
from an operator product expansion 



h(r 1 )...h(r k ) = Y,\ri2\~ kXh+Xc 
o 
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This identity is nothing but a consistency relation for 
the height correlations. Inserted in (^|), it expresses any 
n-point function as a sum of (n — k + l)-point func- 
tions in the limit | r ^ 1 <§C | r^^ | <C £ t (i,j — l,...,k 
and I = k + 1, . . . , n) . The notion of an operator prod- 
uct expansion is familiar in field theory jl^] and has re- 
cently been applied successfully to noncquilibrium sys- 
tems p^5]Jl^| . The sum on the r.h.s. runs over all local 
scaling fields 0(r). Each term contains a dimensionlcss 
scaling function (a simple number for k — 2) and a 
power of |r 12 1 given by the scaling dimensions xq and 
Xh = —x (such that the overall dimension equals that of 
the l.h.s.). The field Ok with the smallest dimension, Xk, 
determines in particular the asymptotic behavior of the 
/c-point functions as £t — > oo, 



<Mri)...fc(r fc )>t~<0fc>t~£t" 
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The amplitudes (Ok)t — (Cfe(r))t diverge with £ t , i.e., 
Xk < |17|. They measure the global roughness, which 
increases as the surface develops higher mountains and 
deeper valleys. Local surface properties should, however, 
behave quite differently. For example, the gradient cor- 
relation functions are assumed to have a finite limit 

lim (Vh(n) . . . Vh(r n )) t = (Vft(ri) . . . Vft(r„)) ■ (10) 



By writing h(vi) — h(r ■) = f 4 ds • V/i(s), the same prop- 
erty follows for the height difference correlation functions 
(n™=i[M r ») — M r 0])*j m particular for the moments (|J). 
This implies a feature familiar from simulations: one can- 
not recognize the value of £t from snapshots of the surface 
in a region much smaller than £ t . 

The operator product expansion (||) induces an expan- 
sion for the gradient field v = V/i of the form 
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with new scaling functions and the dimension x v — 
— X + I- (Both sides of (|ll| ) are tensors of rank k whose 
indices are suppressed.) Hence, the stationarity condi- 
tion (|Io| ) allows only two kinds of terms in (||): 
(a) singular terms involving fields 0(r) with a non- 
negative scaling dimension such as 1 (the identity field), 
(Vft.) 2 (r), etc. Any such field that appears in (|J) will also 
appear in (11) and generate a term (v(ri) . . . v(rfe)) t ~ 



(0) t ~ ^ xo that violates © if x Q < ||. 
(b) regular terms, where the coefficient |r 12 | _ ' £;l: ' 1+:ro 
is a tensor of rank N in the differences rii (i = 2, . . . , k). 
Such terms do not violate ([h]) since they have a vanishing 
coefficient Cf in (0) for N < k. The leading (N = 0) 
term involves the (normal-ordered) field Ofe(r) = h , (r) 
and governs the asymptotic singularity (^J); the higher 
terms correspond to fields with k factors h(r) and N 
powers of V. These composite fields have dimensions 



Xk,N = -kx + N . 
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It is useful to introduce the (normal-ordered) vertex 
fields Z q (r) = exp[qh(r)], which are the generating func- 
tions of the fields h k (r) . Eq. (||) then implies the operator 
product expansion 
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where C° = Cf +i (0, . . . , 0, r 12 /|r 12 |, . . . ,r 12 /|r 12 |) with 
the first k arguments equal to and W4 = c?i |ri 2 | x [p^| . 
Subleading singular terms (with positive-dimensional 
fields O) and regular terms (with fields containing 
height gradients) are omitted. The vertex n-point 
functions {Z qi {r\) . . . Z qn (r n )) t behave asymptotically as 
exp(£* JTJiLi Qi)- If J2iQi = 0; they have a finite limit 
(Z qi (ri) . . , Z-q 1 ...-q n _ 1 {r n )), Since these are precisely 
the vertex correlators that generate the height difference 
correlation functions and since (IK})) is analytic in the g^, 
this leads back to the stationarity condition (|ITJ) . 

The operator product expansion (^) with the linear 
dimensions (fi"^) is at the heart of the field theory for 
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Kardar-Parisi-Zhang systems. It is instructive to com- 
pare this theory with models of turbulence. Burgers' 
equation (||) with force correlations 



rj(r, i)7?(r', t') = a 2 R 2 S(t - t')A(\r - r'\/R) (14) 

over large distances R develops multiscaling: for example, 
the longitudinal velocity difference moments 
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have a fc-dependent singular dependence on 1*12 1 and R 
in the inertial scaling regime a <C |ri2 1 -C R [2(]|§- Sim- 
ilar multiscaling is present in Navier-Stokes turbulence. 
Kolmogorov's famous argument predicts the exact scal- 
ing dimension of the velocity field, x v = —1/3, from 
dimensional analysis ^l}. This determines the scaling 
of the third moment in ( |l5| ) since X3 = 0. The higher 
exponents £4, £5,... < cannot be obtained from di- 
mensional analysis. Assuming the existence of an oper- 
ator product expansion (11), the term ( |l5| ) is generated 
by the lowest-dimensional field Ok with a singular co- 
efficient |2^]. Multiscaling thus implies the existence of 
a (presumably infinite) number of composite fields with 
anomalous negative dimensions. For the velocity vertex 
fields exp[qv(r)] of Burgers turbulence in one dimension, 
Polyakov has conjectured an operator product expansion 
similar to ( |l3| ) p5] . It is no longer analytic in qi , (72 and 
thus consistent with multiscaling. The distinguishing fea- 
ture of Kardar-Parisi-Zhang surfaces is the absence of 
multiscaling p3fl . Notice that the resulting properties 
( |l2| ) and ( |i"3| ) have been derived solely from the assump- 
tions (||) and ([l0|) without using Eq. (|l|) explicitly. 

To establish the consistency of the operator product 
expansion with the underlying dynamic equation, one has 
to construct correlation functions that remain finite in 
the continuum limit a — > 0. With the probability distri- 
bution (0), the height correlations (§) develop singulari- 
ties dictated by their normalization in the linear regime 
(|ry| -C a). The existence of a well-defined asymptotic 
scaling regime for |r,j| ^> a implies that these singulari- 
ties can be absorbed by a change of variables 

h{r) -» Z h (a/r ) h{r) , t - Z t (d/r ) t (16) 

such that the "renormalized" correlations @ satisfy nor- 
malization conditions independently of a at some meso- 
scopic scale tq JlJJl^]. The Z-factors have the asymp- 
totic behavior Z^ ~ (a/ro) x ~ Xo and Z t ~ (a/ro) z ~ z ° as 
d/ro — * 0, where xo = (2 — d)/2 and zq = 2 are the ex- 
ponents in the linear regime. Of course, I do not assume 
perturbative renormalizability (i.e., that the Z-factors are 
analytic functions of Aq). Since the scaling dimensions 
( |l2| ) are linear in k, the renormalization ( |l6| ) also removes 
the singularities from correlations of the fields h k (r) and 
Z q (r) , ensuring a finite limit of the coefficients C in (^), 
Q, @ and of the amplitudes (O k )t in (|). The sub- 
stitution (|l6|) also leads to new coefficients in (j|) and (j^): 
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Galilei invariance is expressed by the asymptotic scale 
invariance of the dimensionlcss coupling rg°A, while the 
other coefficients become irrelevant as a/ro — > 0. How- 
ever, as explained in ref. JL5| for Burgers turbulence, the 
equation of motion for the renormalized correlation func- 
tions is quite subtle due to anomalies dictated by the op- 
erator product expansion. To exhibit the anomalies for 
the height correlations, I introduce the smeared vertex 
fields (r) = exp[g J dr' 8 a (r — r')h(r')] (where S a (r) is 
a normalized function with support in the sphere |r| < a) 
and the abbreviations Zf = Z q . (r^), Z.- L = Z qi (vi). Using 
@, (0), and (|l7|), it is straightforward to derive 



d t {zi . . . zi) t = £ n (zt . . . JZ? . . . Zl) t , (18) 



where JZf = te<r 2 <5 a (0) + J(r,)]Zf . The singularity 
structure of the current is determined by (||) and ( |l7j ) : 

JZf = g*Zi + a 2 *- 2 ^ Oko/^g^j Z t + 0{a x ) (19) 

for a, a — > with a/a kept constant. The field Z q (r) = 
(yh) 2 Z q (v) denotes the finite part of the operator prod- 
uct (Vh) 2 (r)Z£(r) for a -*• 0, and Z, = Z qt {r t ). The 
finite dissipation term (V 2 h)Z qi (ri) becomes irrelevant 
in this limit since v ~ a x . The singular part of ([l9]) is 
a power series in qi with asymptotically constant coeffi- 
cients ci = <r* 2 a d 5 a (0) + v*c± t i +g*C2.i and Ck = v*c\.k + 
g*C2,k for k = 2,3, The terms of order a^ 2 " 1 "*^* -2 orig- 
inate from operator products V 2 /i(ri)/i(r' 1 ) . . . h(r' k ) ~ 1 
and (V/i) 2 (ri)ft,(r' 1 ) . . . h(r' k ) ~ 1; their respective coef- 
ficients ci fe and C2,fc are integrals over the scaling func- 
tions in (g) and the regularizing functions 5 a {^i ~ r j)- Of 
course, divergent terms have to cancel so that Eq. ( jTj ) 
has a finite continuum limit 



d t (Z x ... Z n ) t = 9,: (Zi . . . JZ, ...Z n ) t 
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with j7Zi = lim a ^o i7Zf. For generic values of x, this im- 
plies JZi = g*Zi. However, if x satisfies the condition (||) 
for some integer ko, the dissipation current contributes 
an anomaly: 



JZi = g*Zi + v*ci tko q i °Z i 



(21) 



Eqs. ( p0| ) and (|T]) govern in particular the stationary 
state of the surface. For d = 1, the stationary height 
distribution is known, P ~ exp[— [a 2 /v) J dr(Vh) 2 ]. It 
equals that of the linear theory, thus restoring the up- 
down symmetry h(r) — (h) t — > —h(r) + (h) t broken by 
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the nonlinear term in ([!]). The exponent x — 1/2 satis- 
fies (|J) with fco = 2 but the up-down symmetry forces 
the anomaly to vanish (ci.2 = 0). In higher dimen- 
sions, this symmetry is expected to remain broken in 
the stationary regime. The surface has rounded hilltops 
and steep valleys, just like the upper side of a cumu- 
lus cloud p3]. Hence, the local slope is correlated with 
the relative height, resulting in nonzero odd moments 
((V/i) 2 (i"i)[/i(ri) — h(r2)] k ). However, this is consistent 
with Eqs. (EOh and (Ell) only for odd values of fco, where 



(Z 9 (n)Z_,(r a )> - (Z_ 9 (n)Z g (r 2 )> 

= -(u*/g')c l , ko q k °{Z q {T 1 )Z- q {T 2 )) (22) 

and hence for odd values of k > fco 

((V/O a (ri)[Mri)-M*a)]*> 

= -(v*/g*)c 1 , ko ([h(ri)-h(r 2 )] k - k °). (23) 

The directedness of the stationary growth pattern thus 
requires a nonzero anomaly c± /«„ with an odd integer fc . 
The roughness exponent is then determined by Eq. (|^). 
The values fco = 3 for d = 2 and fco = 5 for d — 3 give the 
exponents quoted above, in reasonable agreement with 
numerical results Jl(],|l3| . 

In summary, the scaling of growing surfaces has been 
determined by requiring consistency of the effective large- 
distance field theory subject to a few phenomenologi- 
cal constraints. The Galilei symmetry of the dynamic 
equation conspires with these constraints to allow only 
discrete values of the roughness exponent in two and 
three dimensions. The underlying solutions of Eq. (|l|) 
are distinguished by a dynamical anomaly in the strong- 
coupling regime: the dissipation term contributes a finite 
part to the effective equation of motion (20) despite being 
formally irrelevant. The anomaly manifests itself in iden- 
tities like ( ^3| ) between stationary correlation functions. 
The quantization rule (||) is analogous to the exact Kol- 
mogorov scaling of the third velocity difference moment 
in Navier-Stokes turbulence. The deeper reason for this 
rigidity is yet to be explained. 

I am grateful to D. Wolf for useful discussions. 
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